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ABSTRACT 

The degree of beaming in Gamma-Ray Bursts (GRBs) is currently unknown. The 
uncertainty in the 7-ray beaming angle, 9^, leaves the total energy release (oc 2 ) and 
the event rate per galaxy (oc 0Z 2 ) unknown to within orders of magnitude. Since the 
delayed radio emission of GRB sources originates from a mildly relativistic shock and 
receives only weak relativistic beaming, the rate of radio-selected transients with no 
GRB counterparts can be used to set an upper limit on 9^ 2 . We find that a VLA 
survey with a sensitivity of ~ 0.1 mJy at 10 GHz could identify ^2x 10 4 (#b/10°) -2 
radio afterglows across the sky if each source is sampled at least twice over a period of 
one month or longer. From the total number of ^ 0.1 mJy sources observed at 8.44 
GHz and the fraction of fading sources at 1.44 GHz, we get the crude limit #b ^ 6°. 

Subject headings: gamma rays: bursts 

1. Introduction 

Gamma-Ray Burst (GRB) sources were discovered historically in 7-rays because they are 
rare and hence require continuous monitoring of large areas of the sky. The necessary all-sky 
monitoring program was first made feasible in the 7-ray regime, hence the name GRBs. By 
now, the prompt 7-ray emission is known to be followed in most cases by delayed emission in 
the A"-ray (Costa et al. 1997), optical (van Paradijs et al. 1997; Bond 1997a,b; Fruchter et al. 
1998a) and radio (Frail et al. 1997a,b; Frail et al. 1998a,b; Kulkarni et al. 1998b; Taylor at al. 
1998) bands. Given the recent discovery of afterglows, it is now timely to explore the feasibility of 
searching for GRB sources at longer wavelengths than traditionally attempted. The importance 
of complementary searches is highlighted by the possibility that some sources might be 7-ray faint 
for geometric or physical reasons. Such sources can only be found at long wavelengths, e.g. in 
radio surveys. 

All past discoveries of GRB afterglows were selected by the detection of 7-rays. If the 7-ray 
emission is beamed, then there should be a population of radio afterglows which are 7-ray faint 
(Rhoads 1997). This follows from the fact that the bulk Lorentz factor of the emitting material 
in GRB sources declines with time and is only of order unity when the late radio emission takes 
place (Waxman, Kulkarni, & Frail 1998), hence making the effect of relativistic beaming weak at 



- 2 - 



that time. In popular synchrotron models (e.g. Meszaros & Rees 1997; Waxman 1997a, 1997b), 
the gradual drift of the peak flux to lower frequencies as a function of time is caused by the decline 
of the Lorentz factor of the emitting material, and so the minimum solid angle into which the 
radiation could be focused increases with time. If the 7-ray emission is beamed into a fraction /b 
of the sky, then the rate of radio afterglows could be larger than the measured rate of GRBs by 
up to ~ f^ 1 . For relativistic expansion with a Lorentz factor T, the radiation observed at a given 
time t originates from a conical section of opening angle ~ 1/T(i) of the fireball. The smallest 
value possible for /b is therefore ~ irT~ 2 /4ir = 1/4T 2 , where T is the Lorentz factor at the time 
of the 7-ray emission. Since typical GRB sources have T ;> 10 2 (e.g., Fenimore, Epstein, & Ho 
1993; Woods & Loeb 1995) while the radio afterglow emission occurs when T ~ 1, the rate of 
radio afterglows could be boosted by ^ 4 order of magnitudes relative to the measured GRB rate. 
The total energy associated with a GRB explosion is proportional to /b and is similarly uncertain. 
Thus, the determination of /b is of crucial importance for identifying the likely astrophysical origin 
of GRBs. 

The evolution of a jet resembles that of a spherical fireball (i.e. the jet behaves as a conical 
section of a spherical fireball), as long as the expansion Lorentz factor is larger than the inverse of 
the jet opening angle, T > Or 1 . The smooth power-law decline of optical afterglows, observed over 
a time scale of days to months for GRB970508 and GRB970228 (Livio et al. 1997; Sokolov et al. 
1998), suggests that these fireballs behaved as if they were spherically symmetric on angular scales 
ranging from 1/T(lday) ~ 1/8 to 1/r (1 month) ~ 1/3 (Rhoads 1997, 1998; Waxman 1997a). 
However, there is still a missing gap in afterglow observations on the time window of minutes to 
several hours following a GRB. During this time the fireball decelerates from T ~ 100 to T ~ 10, 
and so there is uncertainty about the structure of the fireball on angular scales ~ 0.01-0.1 radians. 
It is still possible that the highly relativistic expansion with T ~ 10 2 is restricted only to a very 
small angular diameter, 6*b ~ 1/r = 0.6°. In fact, the popular GRB scenarios of binary coalescence 
of compact stars or failed supernovae favor strong collimation over spherical expansion (Woosley 
1998; Fryer et al. 1998; Meszaros et al. 1998). The search for long- wavelength transients with no 
GRB counterparts can set important constraints on 9^ and hence on these models. 

Current long-wavelength surveys only place weak limits on the beaming factor. In particular, 
existing supernova searches imply that at most ~ 10% of all transient optical events cannot be 
classified as having classical supernovae lightcurves (Kirshner 1998, private communication). Since 
the supernova rate per galaxy is ~ 4 orders of magnitude higher than the GRB rate (Wijers et 
al. 1997), the enhancement factor due to beaming in the optical band must be 10 3 (see Woods 
& Loeb 1998 for a more detailed discussion). The enhancement factor could still be an order of 
magnitude higher in the radio than in the optical. Previous radio experiments such as CLFRT at 
0.151 GHz (Dessenne et al. 1996), and COBE at 90 GHz (Ali et al. 1997) do not provide useful 
limits on this factor, as they achieved only a sensitivity of 10-100 Jy which is well above the 
characteristic mJy radio flux of GRB afterglows. 

In this Letter, we show that a 0.1 mJy radio survey at the VLA [analogous to the FIRST 



survey (Helfand et al. 1996)] can provide strong constraints on 9^. The surveyed sources need to 
be sampled at least twice over a timescale of order a month or longer, since short term variability 
may also be caused by scintillations of steady pointlike sources. 

In our calculations, we assume that the GRB rate is proportional to the cosmic star formation 
rate and neglect the small fraction (<^ 10%) of GRBs which might be related to local radio 
supernovae (Bloom et al. 1998). We model the time-dependent luminosity of the radio afterglows 
based on observational data and allow for a scatter in their peak luminosities. Our model is 
described in detail in §2. The number counts of radio afterglows that it predicts for radio surveys 
are calculated in §3. Finally, §4 summarizes our main conclusions. 

2. Statistics of Radio GRB Afterglows 

The isotropy of GRBs (Meegan et al. 1993; Briggs at al. 1993) and the flattening of their 
number count distribution at faint fluxes suggest that most GRBs occur at cosmological distances. 
This hypothesis has been confirmed by the detection of Fe II and Mg II absorption lines at a 
redshift of z = 0.835 in the optical spectrum of GRB 970508 (Metzger et al. 1997), the inference of 
a redshift z = 3.42 for the host galaxy of GRB 971214 (Kulkarni et al. 1998a), and the detection 
of emission and absorption lines at z = 0.966 for the host of GRB980703 (Djorgovski et al. 1998). 

The energy scale of cosmological GRBs, ~ 10 53 /b ergs (Kulkarni et al. 1998a), corresponds to 
the rest mass energy of a fraction of a solar mass, and implies a link between their energy budget 
and compact stars. Indeed, the most popular models for the origin of GRBs relate them to compact 
stellar remnants, such as neutron stars or black holes (e.g. Paczyhski 1986, 1998; Eichler et al. 
1989; Narayan et al. 1992; Usov 1992; Woosley 1993; Rees 1998; Popham, Woosley, & Fryer 1998; 
Meszaros et al. 1998) . Since these remnants form out of massive stars not long after their birth, 
it is reasonable to assume that the GRB rate traces the star formation rate without a significant 
delay. This scenario was investigated by Wijers et al. (1997), who derived a best-fit constant 
of proportionality between the GRB occurrence rate, R(z), and the star formation rate per 
comoving volume, p(z), based on the requirement that the former fits the observed number count 
distribution of GRBs. The cosmic star formation rate as a function of redshift has been calibrated 
based on observations of the U and U-band luminosity density evolution in the Hubble Deep Field 
(Madau et al. 1996; Madau 1996; Madau 1997; Madau, Pozzetti, & Dickinson 1998). The fit made 
by Wijers et al. (1997) yields a local GRB rate of R(z = 0) = (0.14 ± 0.02) x 10~ 9 Mpc" 3 yr- 1 in 
an f2 = 1, A = 0, Hq = 70 km s^ 1 Mpc -1 cosmology. 

We model the frequency and time dependence of the afterglow luminosity, L u (t), using the 
simplest unbeamed synchrotron model (e.g. Waxman 1997a, 1997b), with a luminosity per unit 
emitted frequency v of 
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where v m 

(t) = 8.8 x 10 2 (1 + z) 1 / 2 (t/month)- 3 / 2 GHz, t is time at the source frame, and z is 
the source redshift. The observed frequency v is related to the emitted frequency v through 
v = P/(l + z). The spectral index a is chosen to have the values a\ = 1/3 for v < u m and = 0.7 
for i/ > i/ m , so as to match the temporal decay slope observed for GRB 970228 (Galama et al. 
1997; Fruchter et al. 1998b) and GRB 970508 (Galama et al. 1998). 

We consider a population of afterglow sources characterized by a total comoving rate per 
unit volume, R(z), and by a peak flux F Um {z,L Um ) = L Vm {l + z)/AirDl(z), where D^(z) is the 
cosmology-dependent luminosity distance. For consistency with the numbers derived by Wijers et 
al. (1997) we assume Q = 1, A = 0, h = 0.7, and D L (z) = (2c/H )(l + z- V^Tz). Note that if 
the radio flux of GRB sources is strictly proportional to their 7-ray flux, then our number count 
results are independent of the choice of cosmological parameters to zeroth order; this follows from 
the fact that Wijers et al. (1997) calibrate their results based on GRB number count data. 

We allow for a spread in the peak luminosity L Um of the afterglows by using a log-normal 
probability distribution (and thus minimizing the number of free parameters), 

P(L )dL - 1 cxpf N^ m )-1°(^)] 2 ]^ (2) 

The mean of this distribution equals (L Um ) = L*exp (<t 2 /2). We normalize its value by matching 
the mean peak flux, 

lF f °° dL Um P(L Um ) J °° dz[R(z)/(l + z)]F Um (z, L Um )dV c 

{ Vm) " fo°° dL Vm P{L Vm ) / °° dz[R(z)/(l + z)]dV c ' 1 > 

to its observed value in radio afterglows. The observed 8.44 GHz lightcurves of GRB970508, 
GRB980329, and GRB980703 yield (F Um ) w 0.7 mJy. For the above cosmology, this implies 
{L U J w 1.5 x 10 31 /if 2 erg sec" 1 Hz" 1 . 

The calculation of the afterglow number counts should also incorporate the fact that not 
all GRBs have afterglows at radio frequencies. The fraction / ra( jio of radio-loud GRBs is still 
uncertain, but the observation that 4 out of ~ 15-17 well-localized GRBs (namely, GRB970508, 
GRB980329, GRB980519 and GRB980703) were detected in the radio (D. Frail 1998, private 
communication), implies that / ra( jio — 25%. Also, the afterglow radio flux might have a break in 
its power-law evolution and decay more rapidly than predicted by equation (|l]) after a particular 
time (e.g., during the non-relativistic phase of the shock hydrodynamics). We artificially introduce 
a free parameter, t C utoffj after which we set the emitted radio flux to zero. Current data does not 
allow an empirical calibration of Cutoff an d so we leave it as a free parameter. We assume that 
icutoff ^ 3 months since a simple power-law decline was observed in GRB970508 over 3 months 
(Waxman et al. 1998). 

The late phase of the radio emission, several months after a GRB event, originates from 
semi-relativistic material with T ~ 1 and receives only weak relativistic beaming. However, the 
7-ray emission is emitted by material with T ^ 10 2 and could be confined to a cone with an 
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angular diameter #b *C 1. In this case, the rate of radio afterglows would be enhanced by a factor 
/- 1 = 47r/7r(0 b /2) 2 = 5.25 x 10 2 (<9 b /10°)- 2 , relative to the unbeamed case. 

Typical radio surveys do not provide continuous monitoring of the sky, but rather a sequence 
of "snapshots" of each resolution element on the sky. In observations where the minimum 
detectable flux of a source is F u , the total number of events across the sky that are brighter than 
F u at observed frequency v is given by 

r-oo 

N(> F„; v) = / b - Vradio / dL Um P{L Um )N'{> F v \ u, L Um ) (4) 

Jo 

where 

rzi im (L Vm ) fjy 
N'{> F v ;v,L Vm ) = J R(z)t*(z,F u ,v,L Um )^dz . (5) 

Here dV c = 4(c/# ) 3 (l + z- y/T+z) 2 (l + z)~ 7 / 2 dndz is the comoving volume element within 
a solid angle dd and redshift interval dz, and z\\ ni (L Um ) is found from the algebraic relation 
Fv m {z\\ m -, L Vm ) = F v . The time t*(z,F v ,i/) represents the duration over which an event at a 
redshift z is brighter than the limiting flux F v at the observed frequency v, 



(l+z)u\- 2/?> 



U(z, F„, is, L Um ) = min <^ t cuto ff, 



where u mo = 8.8 x 10 2 (1 + z) 1 / 2 GHz month 3 / 2 . 



2_ 2_ 

F v \ 3a 2 / F v \ 3a 1 



F u m (z,L Um )J \F Um (z,L u , t 
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3. Constraints from Radio Surveys 

The VLA FIRST survey monitored 1550 degrees of sky with a sensitivity of 1 mJy at a 
frequency of 1.5 GHz (Helfand et al. 1996). The intervals between the observations of each source 
ranged between 3 minutes and 3 weeks. Unfortunately, the properties of this survey were not 
optimal for the purpose of identifying GRB afterglows. Here we consider a hypothetical survey 
which is optimized for this purpose. Such a survey would use a higher frequency in order to 
avoid the suppression of the afterglow flux by synchrotron self-absorption at the source (typically 
observed at <^ 5[(1 + z)/2]~ 1 GHz). Moreover, since short-term variability could also be caused by 
scintillations, the surveyed sources should be monitored at least twice over an extended period of 
order a month or more. After a few months all the radio-loud afterglows would fade significantly 
and thus be easily distinguishable from other variable source populations. The variability level of 
~ 0.1-1 mJy sources can be determined accurately, given the high sensitivity of the VLA. 

Figure 1 shows the number of radio afterglows that are expected to be found in an all-sky 
survey at 10 GHz, as a function of its flux threshold. The plotted number depends on a 
normalization factor involving the 7-ray beaming angle #b in units of 10°. We assume / ra dio = 25% 
and show results for different values of a while keeping t C utoff fixed at 6 months. Figure 2 shows 
the same for a — > but for different values of t cuto g. Typically, a 0.1 mJy survey similar to FIRST 
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but at 10 GHz would identify ^2x 10 4 (#b/10°) 2 afterglows across the sky if the 7-ray emission 
from GRBs is collimated to within an opening angle b - 

Lower limits on the value of 9b can already be placed by existing radio surveys. The total 
number of radio sources observed on the sky provides a robust upper limit on the number of GRB 
afterglows. This limit can be improved by considering only those sources which are unresolved 
and at the same time variable over a time scale of a month or longer. Still tighter constraints can 
be obtained by considering only the subset of all unidentified sources that after a sufficiently long 
time, ;> 6 months, fade away below the detection threshold. From a deep VLA survey, Windhorst 
et al. (1993) derived a maximum-likelyhood fit to the source number counts at 8.44 GHz in the 
flux range of 14.5//Jy-1.5mJy. By comparing a Westerbork and a VLA survey of the same field 
at 1.4 GHz, Oort & Windhorst (1985) found that about 10% of the sources were variable on a 
time scale of a year, and only ~ 3% were detected in one survey but not in the other. Thus, only 
~ 3% of the sources in this survey qualify as candidates for GRB afterglows. The sub-mJy sources 
were typically found to be unresolved. Better angular resolution could in principle limit further 
the number of afterglow candidates by resolving other source populations. 

Figure 3 shows the ratio between the predicted number of GRB afterglows N(> F v ) and the 
total observed number of radio sources N b s (Windhorst et al. 1993) at 8.44 GHz. We denote the 
fraction of variable sources which fade away after a year by /fade; and the fraction of pointlike 
(unresolved) sources by / p0 mt- Based on the observed abundance of unresolved variable sub-mJy 
sources at 1.4 GHz (Oort & Windhorst 1985), we set /f a( je x /point = 3%. The ratio between 
the predicted afterglow counts and the observed radio counts peaks around 0.4-0.6 mJy. From 
the requirement that this ratio be smaller than unity, we get the lower limit 9b ^ 6°. This limit 
is based on the variability statistics of only ~ 10 2 sources (Oort & Windhorst 1985) and hence 

— 1/2 

suffers from large statistical uncertainties (min{^b} oc / fad g ). Surveys with a larger coverage of 
the sky are necessary in order to firm up the statistical significance of this limit. 

4. Conclusions 

We have found that a 0.1 mJy all-sky survey at 10 GHz should identify ~ 2 x 10 4 (6> b /10 o )~ 2 
GRB afterglows, and could therefore place important constraints on the GRB beaming angle, 
6> D . In fact, the total number of ^ 0.1 mJy sources on the sky at 8.44 GHz (Windhorst et al. 
1993) and the fraction of fading unresolved sources at 1.4 GHz (Oort & Windhorst 1985), already 
yield the crude lower limit 9b ^ 6°. This result is only weakly sensitive to the width of the 
afterglow luminosity function (see Fig. 1), and is mainly uncertain due to the limited size of the 
sample of sub-mJy sources that had been monitored for long-term variability (Oort & Windhorst 
1985). Nevertheless, the derived constraint is interesting, given the fact that the 7-ray emission 
in GRB sources originates from material with a bulk Lorentz factor V ^ 10 2 and could have been 
collimated to within an angle b ~ 1/r <; 0.6°, narrower by more than an order of magnitude than 
our limit. In comparison, radio jets in quasars are collimated to within ~ 15° and possess Lorentz 
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factors r ^ 10 (Begelman, Blandford, & Rees 1984). 

A future VLA survey, optimized to search for radio afterglows by monitoring sub-mJy source 
variability at ~ 10 GHz over a timescale of several weeks to several months, could improve the 
above upper limit on 9^ considerably. Alternatively, such a search might identify a new class of 
radio afterglow events which fade away after a few months and have no GRB counterparts. 

We thank Dale Frail for useful discussions. This work was supported in part by the NASA 
grant NAG5-7039. 
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Fig. 1. — Number of radio-loud afterglows that would be detected in an all-sky radio survey as a 
function of its flux detection threshold F v at 10 GHz. The different lines correspond to different 
values of the width of the GRB luminosity function, a [see Eq. (2)]. The fraction of the radio-loud 
afterglows is assumed to be / ra dio = 25% and the afterglow emission is truncated after a time 
^cutoff = 6 months in all sources. The normalization factor of the number counts involves the 
opening angle of the 7-ray emission, 6>b, in units of 10°. 
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Fig. 2.- 



Same as in Figure 1, but with a fixed cr 



and different values of t C utoff • 
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Fig. 3. — The ratio between the predicted cumulative number of GRB afterglows iV(> F v ) and 
the observed cumulative number of radio sources iV bs (Windhorst et al. 1993), as a function of 
the flux threshold F u at 8.44 GHz. The observed number is uncertain by ±40%. The fraction of 
unresolved fading sources in the observed population is assumed to be /fade/point = 3%, based on 
1.4 GHz data (Oort & Windhorst 1985). Model parameters are chosen as in Figure 1, with a — > 0. 
The condition that the plotted ratio be less than unity at any flux threshold sets a lower limit on 
the 7-ray beaming angle 6^ ;> 6°. 



